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We consider spin-half fermionic atoms with isotropic Rashba spin-orbit coupling in three direc- 
tions. The interatomic potential is modeled by a square well potential. We derive the analytic form 
of the asymptotic wave-functions at short range of two fermions in the subspace of zero net mo- 
mentum and zero total angular momentum. We show that the spin-orbit coupling has perturbative 
effects on the short range asymptotic behavior of the wave-functions away from resonances. We 
argue that our conclusion should hold generally. 



I. INTRODUCTION 

The study of dilute unitary atomic Fermi gases has sig- 
nificantly and substantially extended our knowledge on 
strongly interacting many-body systems [TJ. Theoretic 
advances in this direction have been propelled by ob- 
serving the short-range asymptotic behavior of the wave- 
function of the systems [2H1]- In such ultra-cold dilute 
gases, when the separation r between two fermions of 
different species is much smaller than the mean inter- 
particle spacing d, but bigger than rrj, the range of the 
interatomic potential, the pair wave-function of the two 
fermions is mainly s-wave and has the asymptotic form 
ip s (r) = 1 — a s /r; when a similar situation happens to 
two fermions of same species, the pair wave-function is 
mainly p-wave and its radial part has the asymptotic 
form ip p (r) — k(r — u p /r 2 ). Here k is the relative wave 
vector between the two fermions and of order 1/d, a s is 
the s-wave scattering length and u p is the p-wave scat- 
tering volume. Thus, when considering the short-range 
correlations, the p-wave part can be safely neglected com- 
pared to the s-wave part unless u p becomes divergent. 
Based on this observation, a collection of remarkable re- 
lations regarding the short-range correlations and various 
physical observables have been derived [2H1]- 

On the other hand, recent successful engineering of 
synthetic gauge field adds another important dimension 
to atomic gases [5HTU] . Raman processes couple different 
hyper-fine states of atoms and realizes a model with an 
effective spin-orbit coupling for spin-half particles when 
the hyper-fine states other than the lowest two are adi- 
abatically eliminated. This success raises the prospect 
of using ultra-cold atomic systems to study and simulate 
spin-orbit coupling physics [IT], which has excited great 
interest in condensed matter physics [TH [T3] . 

The introduction of the spin-orbit coupling to ultra- 
cold atomic gases also raises another interesting ques- 
tion: How would the short-range asymptotic behavior 
of the wave-function of atomic gases be modified? At- 
tempts to answer this question include: Cui studied the 
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two-body problem of fermions with symmetric Rashba 
spin-orbit coupling in three directions and focused on 
properties of the bound states. The discovery that the 
usual s-wave pseudo-potential is still applicable in cer- 
tain regimes indicates that the effects of the spin-orbit 
coupling on the short-range asymptotic behavior of the 
wave-function can be perturbative [14]. Later Ref. [15] 
presented a general argument that since the spin-orbit 
coupling strength, as usually realized in the experiments, 
is much smaller than the energy scales associated with in- 
teratomic potentials, the wave-function inside the range 
of the interatomic potential should remain intact to zero 
order. The asymptotic behavior of the wave-function out- 
side the interatomic potential range which is determined 
by the one inside the range should stay unchanged as 
well. Re fence [16) reached a similar conclusion by an 
argument in which unitary transformations are imple- 
mented to eliminate the spin-orbit coupling from the ki- 
netic part of the Hamiltonian. However, besides the case 
with Rashba spin-orbit coupling in one direction [16] , the 
general arguments mentioned above await justification 
from explicitly worked out examples. 

In this paper, we study the problem of two spin-half 
fermions with isotropic spin-orbit coupling in three di- 
rections. The single-particle Hamiltonian of the fermions 
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with m the mass of the atoms. Without loss of generality, 
we assume the spin-orbit coupling strength to be positive, 
A > 0. We model the interatomic potential by an attrac- 
tive square well potential, V(r) = —Vo9(ro — r), where tq 
is the potential range and Vq > 0. As of experimental in- 
terest, we assume Aro -C 1. We calculate the scattering 
and bound state wave-functions of the two interacting 
fermions in the subspace of zero net momentum and zero 
total angular momentum (sum of the spin and orbital 
angular momentum). We derive analytically the asymp- 
totic behavior of the wave-functions outside the range 
of the square well potential. Away from resonances, the 
modification of the asymptotic behavior due to the spin- 
orbit coupling is shown to be perturbative. We explain 
this perturbative effect by inspecting the wave-functions 



2 



inside the potential range and argue that our conclusion 
is valid for generic situations with spin-orbit coupling. 



II. THE SCHRODINGER EQUATION 



III. SCATTERING STATES 



It is straightforward to show that for the scattering 
states, the wavefunctions are 



The Hamiltonian of two interacting fermions can be 
cast into the form 
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Here K is the net momentum of the two fermions and k is 
the relative one, and r are the relative coordinates. The 
subscript of cr labels for the «th-fermion. From Eqs. ([3| 
and Q , the motion of the center of mass and the relative 
one are coupled together via the spin operators. In the 
following, we focus on the subspace of K = and J = 
L + S = 0, with L the orbital angular momentum and S 
the total spin, where analytic results can be derived. 

The relative wave- function ^(r) in the subspace K = 
satisfies the Schrodinger equation 
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when represented in the spin basis — ItVv^ttill 
,(U + lt)/V2)} for the two spin-half fermions. Since 
each elements of M\ are proportional to the spherical 
harmonic functions Yi jTO (^k) or their complex conju- 
gates, with further contraint J = we have 
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for r < ro, and 
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for r > r . Here ji are the spherical Bessel functions and 
are the spherical Hankel functions, and 

k\ = k - A, k 2 = k + A, (11) 
9l = A - q, q 2 = A + q, (12) 
k = VmE, q ee ^m(E + V ). (13) 



Since there are six coefficients (^4 to F) and four con- 
nection conditions at r = ro for the wavefunctions and 
an overall normalization factor, the wavefunctions can 
be determined up to a coefficient; there are two lin- 
early independent solutions for each energy E. Let 
us choose the ith linearly independent solution $^ as 
$i = [3o{qi~r),ji(qir)} T for r < r . By requiring 
and its first derivative continuous at r = ro, we find the 
corresponding coefficients for $?" 
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The new spinor wave-function $(r) = [V>o( r ). 4 } i{ r )] T sat- 
isfies 



\^(r 2 ~)+>? + m[V(r)-E] 
r z ar \ ar 



-2\{d/dr + 2/r) 
2A (d/dr) 2/r 2 



$( r ) = 0. (8) 



with ki = kiro and <ji = qiTQ. 

The physical meaning of the above coefficients can be 
understood in the following way. Without interaction, 
V(r) = 0, cither [jo(kir),j 1 (k 1 r)} T or [j (k 2 r), ji{k 2 r)) T 
is the free particle solution to Eq. |8]). The interaction 
V(r) realizes mutual scattering between the two waves. 
We can construct two new linearly independent solutions 
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where Sij are the elements of the matrix 



E 2 C\ — E\C 2 
^(E^-E*^ 



^{C* 2 C X -CIC 2 ) 
C\E 2 — C 2 E\ 



(20) 



Since (x) 



jx and ft- 2 ^(a;) ~ e ,a, /x when a; 



00, the wave- function $^ ) describes the process that 
the incoming wave of wave vector k\{k 2 ) is scattered by 
V(r) into the outgoing waves of wave vectors k\ and k 2 . 
The matrix 5 formed by the coefficients is the scat- 
tering S'-matrix. It is straightforward to show that 5 
satisfies the unitary condition S^S = 1. The appearance 
of the ratios between k\ and k 2 in Eq. ( 20 ) is because the 



two wave vectors are of different magnitude. 

We can further obtain the wave-functions such that the 
S-matrix is diagonalized: 
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where [vf,vf] are the eigenvectors of S given by 

Eq. (20) and e 2lS± are the corresponding eigenvalues. 
From Eqs. (14) and (15), one can prove k^E^Ei — 
E{E 2 ) =kl{U£Ci- CfC 2 ). Since the S-matrix can be 
expanded as S = Sol+S x a x +S z a z and the ratio between 
the coefficients S x and S z is real, we can choose the eigen- 
vectors \vf, v 2 ] T to be real. The solutions <I>± are stand- 
ing waves since [$±]* is proportional to $± apart from 
a phase. This complies with the expectation that solu- 
tions to the one-dimensional differential equation ^ can 
be chosen to be real. Equation (21) corresponds to the 
ansatz used in Ref. P3] [cf. Eqs. (31) and (32) therein]. 
The two new phase shifts 5± are the counterparts of the 
s-wave and p-wave phase shifts 5 S and 5 P in the absence 
of the spin-orbit coupling. 

We plot S_(kr ) in Fig. Q and 5+(kr ) in Fig. ^ 
for A = 0.01 and rj = 1 with A = Xtq and 77 = 
\J mVorQ. Generically we find exp[2i(5_(0)] = —1 and 

exp[2z<5+(0)] = 1. The unitarity of (5 (0) is a reminiscent 

of one-dimensional scattering [13]. When both <5_ and 
(5+ are defined in the domain [— ir/2, ir/2], to zero order 
of A, <5_(0) jumps by it where a s switches between — 00 to 
+00, and the slope of 5-(k) in the small k limit changes 
sign where u p flips between — 00 to +00. To analyze the 
properties of 8±, we expand the S'-matrix to first order 
of k 
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where k = kr$. The elements Mij of the matrix M are 
listed in the Appendix. Correspondingly in such limit, 
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+ 2Acos(2r ? )[-2j 7 2 A + A 3 -r ? 2 sin(2A)]}/{2?7 2 A 2 [cos(277) -cos(2A)] +2t]A 2 (j7 2 - A 2 ) sin(2?7)}, (25) 

5' + = - *(Af n + M 22 + M 12 + M 21 ) 

={r/[r/ 2 + A 2 - 2j7 2 A 2 + 2A 4 - (r] 2 + A 2 ) cos(2A)] sin(2?7) + 2A cos(2r/)[A 3 - 2?7 2 A + rj 2 sin(2A)] 

- 2A[A 3 - r/ 2 A - ?7 2 Acos(2A) + r/ 2 sin(2A)]}/{2,] 2 A 2 [cos(2,7) - cos(2A)] + 2t 1 \ 2 (t 1 2 - A 2 ) sin(2?7)}, (26) 
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Sv =(M n - M 22 )/A 

={2?7A[cos(27?) cos(2A) - 1] + {rj 2 + A 2 ) sin(2r/) sin(2A)}/{2?7A 2 [cos(2r;) - cos(2A)] + 2A 2 (r? 2 - A 2 ) sin(2?7)}. (27) 
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FIG. 1: Phase shift 5-(kr ) vs kr for A = 0.01 and r] = l. 



kr„ 



FIG. 2: Phase shift 5+{kr ) vs kr for A = 0.01 and rj = 1. 

The three quantities, 5'_, 5', and 5v, diverge when 

71 ^ [cos(2ry) - cos(2A)] + (rj + A) sin^ry) = 0, (28) 
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where a bound state emerges [cf. Eq. (38)]. This diver- 
gence reflects that bound state energies are the poles of 
the S'-matrix. Away from the singularities, these quanti- 
ties have the asymptotic expansions in A as 
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For our square well potential model, the s-wave scat- 
tering length a s is given by r Q /a s — r)/(r) — tan 77), and 
the p-wave scattering v olum e u p by Up/r^ = 1 — 3(1 — 
rj cot rj)/r] 2 . Equations (|29|) and (30) explain the jumps 



of <5_(0) and the sign change of the slope of 5+ (k) for 
small k. 

Exceptions occur at resonances: at the first set of res- 
onances given by Eq. (41), we find S-(Q) = — <5+(0) « 
—A/2, v = [0, 1] T , and v + = [1, 0] T ; at the second set, 
5_(0) = -S + (0) « 2A 3 /3, v- = [0, 1] T , and v+ = [1, 0] T . 



The asymptotic form of the scattering wave-functions 
in the region r > rg is 
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Note that an overall factor 2e s± sin S± has been omit- 
ted. In the low energy limit k — > 0, when away from 
resonances, using Eqs. (29) to (31), we find 
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At the first set of resonances given by Eq. (|4l|) 
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Note that Eqs. Q to Q satisfy Eq. ([8) with £ = to 
the leading order of A. 



IV. BOUND STATES 



The bound state solutions to Eq. ^ with energy E 
are Eq. Q for r < r and 
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for r > r . Here fc' x = A + in, k' 2 = A — and k = 
y/—mE. The requirement that the wave-function and its 
first derivative are continuous at r = ro gives rise to the 
equation determining the energy E: 

=(A 2 + k 2 + k){q 2 - K 2 )bo(9i)ji(g 2 ) - h(h)h{qx)\ 

+ 2AK0o(<?2)jl(<Zl) + JO («l)jl(?2)] 

-2^[(K+l) 2 + A 2 ]j (gi)jo(g 2 ) 
-2^(A 2 + K 2 )j 1 (9i)jife), (38) 
with k = ktq and (7 = grp. Assuming D' ~ 1, we have 
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In the weak attraction limit rj — > 0, we find 
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there is always a bound state no matter how weak the 
attraction is. This bound state disappears in the absence 
of the spin-orbit coupling since from Eq. ( 40 1 k — > for 
A 0. The energies of bound states other than in the 
weak attraction limit have been studied in Ref. [14] . 

According to Eq. ( 38 ) , a new bound state emerges at 
threshold E = when Eq. (28) is satisfied. In the limit 
A -> 0, Eq. (|28j) reduces to j (rj) = or ,71(77) = 0. In the 
absence of the spin-orbit coupling, the former condition 
is where a bound state emerges in the p-wave channel, 
while the latter is where the s-wave scattering length a s 
is zero. The introduction of the spin-orbit coupling has 
the nonperturbative effect of shifting the onset of a series 
of bound states from where l/a s = to where a s = 0. 
For small A, the critical values of the depth of square well 
potential V c where a bound state forms are given by 
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with m — 1,2,3,... and nonzero £ n (n — 1,2,3,...) 
being the solutions of ji(l n ) = 0. The change of V c due 
to small A agrees with the numerical results obtained in 
Ref. [H]. 

To obtain the asymptotic behavior of the wave- 
function of the most shallow bound state for r > r when 
Vq is close to a critical value V c , let us first look at the 
limit 77 — > B£_= where there is only one bound state. 
From Eq. (40 1, in the limit A — > 0, for the bound state, 
we have R = ry 2 A 2 /3. Substituting this into Eq. (391 and 



expanding in the order 0<K<g;77<§;A<?Cl, we obtain 

= 2arctaaQ [l + C(A 2 )] ~ 2^, (42) 

and thus for r > ro the wave-function of the bound state 
asymptotically is 
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which coincides with the 77 — >■ limit of Eq. ( 33 ) 



When Vq approaches one of the first set of critical val- 
ues of Eq. (Ell), 
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which agrees with Eq. ( 35 ) to the corresponding order 
of r. For Vq close to the second set of critical values of 
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which is the same as Eq. (|36j) to the corresponding order 
of r. 



V. DISCUSSION 



Equations ( 33 1 and ( 34 ) indicate that away from res- 



onances the modification of the asymptotic forms of the 
wave-functions at short range due to the spin-orbit cou- 
pling is perturbative. To see this point clearly, we multi- 
ply Eq. ( 34 1 by A and have 
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apart from certain factors, to lowest order of A, 
have the same asymptotic form as the s-wave- function 
ip s and $^ as the p-wave-function tp p . The perturbative 
modification can be understood by inspecting the wave- 
functions inside the range of the potential. For r < ro, 
in the limit k —¥ and expanded to the first order of A, 
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Thus to zero order of A, $+(r) and $^(r) are the same 
as the s- and p-waves inside the square well potential in 
the absence of the spin-orbit coupling respectively. Since 
the coefficients of the terms linear in A are well-behaved 
away from the resonances, nonzero A adds corrections 
perturbatively. Since $ > (r) and $ < (r) need to connect 
at r — r , the change of the asymptotic form of $ > (r) is 
bound to be perturbative in A as well. 

The above analysis is essentailly the same as the pre- 
vious argument based on the comparison of energy scales 
given in Ref. |15) : Within the potential range, the spin- 
orbit coupling strength is much smaller than 1/Vo and 
\/mVo and therefore generally has perturbative effects 
on the wave-functions in this regime, and so does on the 
wave-functions in the regime r > tq. According to this 
general consideration, we expect the same conclusion ap- 
plies to the cases with nonzero pair net momentum K and 
J ^ 0, and to the ones with generic forms of spin-orbit 
coupling, given additionally K <C 1/ro, y/mVo- Since 
Ref. [15] assumes implicitly u p — 0, the result derived 
there agrees with Eq. (48) while there is no match for 
Eq. {49J. 

Nevertheless, Eqs. (50) and ( p3lj ) signal that the per- 
turbation expansions in A break down when a s — > 
(77 — tan 77 = 0) and l/u p — > (tan 77 = for 77 ^ 0), which 
are the conditions for the emergence of a new bound state 
in the limit A — > 0. At the first set of resonances, the 



perturbatively. At the second set of resonances, extra 



powers of A in the denomenators of ( 36 ) become singular 
as A — > there. This coincides with the perturbative 
and nonperturbative shifts of the onset of bound states 
compared to the cases without the spin-orbit coupling 
mentioned in Sec. IIVI 

Our calculations suggest that for generic spin-orbit 
coupling, no matter for bosons or fermions, the correc- 
tions to the asymptotic form of the wave-functions at 
short range should be perturbative except for at fine 
tuned points where resonances occur. Consequently the 
relations derived based on the asymptotic form [5H3] 
should stand valid correspondingly. 
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Appendix A: Elements of the M matrix 



asymptotic forms ( 35 1 suggests that A can still be treated 



We list the elements of the M matrix here: 

_ e- 2 ^-^r 1 (r ] + A) 2 - e^+^fa - A) 2 + 4t\[ V 2 (i + A) - A 3 ] 

11 ~ 2t7A 2 {77[cos(2t7) - cos(2A)] + [rf- - A 2 ) sin(2r/)} ~' 
_ .2A 2 (2r/ 2 - A 2 )cos(2r/) - 2?7 2 A 2 cos(2A) + 

12 ~ ' 77A 2 {77[cos(2t7) - cos(2A)] + (r? 2 - A 2 ) sin^)} ~ 

M21 =M 12 , (A3) 
M 22 = -M* 1 . (A4) 
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